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We consider certain examples of applications of the general methods, based on geometry and in- 
tegrability of matrix models, described in £Q. In particular, the nonlinear differential equations, 
satisfied by quasiclassical tau-functions are investigated. We also discuss a similar quasiclassical 
geometric picture, arising in the context of multidimensional supersymmetric gauge theories and the 
, AdS/CFT correspondence. 

\D 

O - 

1 Introduction 



I In the first part of this paper [*Q we have discussed the properties of the simplest gauge theories - the matrix 
integrals 

r i r N 

Z = j d*e-*^W = JL j JJ (^ e -t^)) A 2 (^>) (1.1) 



X 



1=1 

and 

N 



Z = J d$d$t cxp (-Iv^,^ = i Jf[ (d 2 Zl e--* v ^) |A(z)f 



(1.2) 



or the so called one-matrix and two-matrix models, where by 

V(<P,<f> r ) = - W(®;t) - W(&;t) (1.3) 

and 



k 

fe>0 



(1.4) 



the corresponding gauge-invariant single-trace potentials (the last one is usually some generic polynomial) are 
denoted, and 



A(0) = n(^-^i) = det||^- 1 



(1.5) 

i<i 

* Based on lectures presented at several schools on mathematical physics and the talks at "Complex geometry and string theory" 
and the Polivanov memorial seminar. 
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stays for the Van-der-Monde determinant. We have reminded in 1 that matrix models Ijl.lfl . (|1.2f> are effectively 
described in terms of Toda integrable systems, and demonstrated, that in planar limit N — > oo of matrices of 
infinite size (or h — > with hN = to = fixed) these integrable systems have nice geometric origin. 

The planar limit corresponds to extracting the leading contribution to the free energy of a gauge theory (in 
particular, of 11.211 ) 

-Zgaugc = exp (String) (1.6) 

which is a certain string partition function (see, e.g. 0), dual to a gauge theory. In particular case of the 
matrix ensembles l|l.l[l or l|1.2fl . within their 1/iV-expansion |3] 

oo 

^string = Yl N 2 - 29 F g (t, t Q ) (1.7) 
9=0 

the free energies Fq = T of the planar matrix models can be identified with the quasiclassical tau- functions 0j , 
or prepotentials of one-dimensional complex manifolds S. These one dimensional complex manifolds or complex 
curves have the form 

In 

y 2 = W'(x) 2 +4f(x) = Hix-Xj) (1.8) 
j'=i 

for the one-matrix model (|l.l|l . and 

z n ~z n + a n z n+1 + a n z n+1 + hi zt ~ zJ = 

l,j£{N.P.) + 



(1.9) 



in the two-matrix case ljl.2JI correspondingly; in the last formula the sum is taken over all integer points strictly 
inside the corresponding Newton polygon, see details in 1 . To complete the geometric formulation, these curves 
should be endowed with the meromorphic generating one-forms dS 

— ydx, one — matrix model 

J (1.10) 

zdz, two — matrix model 

2™ 

and the corresponding prepotentials are defined in the following way: 

S = I dS, 

QT (LID 

ds- 2m i ds 

i.e. as functions of half of the periods of the generating one-form (|1.10|) so that their gradients are given by 
dual periods w.r.t. the intersection form A a o Bp — 5 a/ 3. Integrability of Ijl.llfl is guaranteed by symmetricity 

i£k = 2mTap (L12) 

of the period matrix of S, which is a particular case of Riemann bilinear relations. 

In the case of matrix models and their curves (|1 . S|> . (|1.9|l the period variables (|l.ll|l have meaning of the 
fractions of eigenvalues, located at particular support. The total number of eigenvalues 

to = resp.dS = respdS 

ziri 2tti 

d _F r? + (1-13) 

dt 



= Am I dS 
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and the parameters of potentials l|1.4fl . (|1.3fl 



t k = 7T-T res Po£ k dS, k > 

dT 1 (1-14) 
— = — res Po £ k dS, k > 

are associated with the generalized periods or residues of generating one-forms; in 11.14fl £ is an inverse local 
co-ordinate at the marked point Pq: £(Pq) = oo. Analogously to l|1.12f> the consistency condition for H1.14fl is 
ensured by the following formula for the second derivatives 

while the third derivatives arc given by the residue formula 

d^T 1 f dH I dH J dH K \ 1 ^ ( fa 



dT ld T jd T K 2m ieSd *=° { 'dxdy " ) = 2^ "** V^W^ ) ^ 
where the set of one-forms {dHi} = {dtv a , dQo, dfl^} corresponds to the set of parameters {T/} = {S a ,to,tk} 

by 

ddS , . . 

-— =<iw a , a = l,...,g (1.17) 

oS a 

together with 

^ = dn (1.18) 

dt 

and 

^ = c/n fe , *>1 (1.19) 
In the second line of (jl . lf>|> - for convenience in what follows, we have introduced the meromorphic functions 

, , dHj 
ay 

As functions of their parameters, being (extended) moduli of complex manifolds (|1.8|) . (|1.9|) the quasiclassical 
tau function or prepotential is defined only locally, but in the way consistent with the duality transformations: 
A <-> B 1 in the sense of the corresponding Legendre transform 

T <-> T + ^2s a n a ( 12 i) 

a 

For different choices of the basis cycles the functions T differ merely by notation, and we will not distinguish 
between different functions below when studying equations, satisfied by quasiclassical tau-functions of matrix 
models. 

In the next section we start with studying nonlinear differential equations, satisfied by prepotentials or 
quasiclassical tau functions, which contain the dispersionless Hirota equations and their higher genus analogs and 
the associativity or Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equations. Then we turn to various physical 
examples of the gauge/string duality Ijl.tifl where the considered in methods of "geometric integrability" can 
be applied. We start with the well-known topic of duality between matrix models and two-dimensional gravity, 
formulated along these lines in terms of dispersionless hierarchy and then turn to the geometric picture of the 
AdS/CFT correspondence. 
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2 Dispersionless Hirota equations 



In this section we consider the differential equations, satisfied by quasiclassical tau-function. We first discuss 
in details the rational case, where E is Riemann sphere P 1 and it is possible to write easily the explicit form 
of differential equations for the quasiclassical tau-function, appearing to be the dispersionless limit of the 
Kadomtsev-Petviashvili (KP) and Toda hierarchy Hirota equations. Finally in this section, we discuss the 
analogs of dispersionless Hirota equations for the curve E of higher genus and the WDVV equations. 

2.1 Geometry of dispersionless hierarchies 

Consider a generating function: the Abelian integral S = J P ydx, corresponding to the generating form (|1.10|) 
of the one-matrix model. From the formulas (|1.14|) one finds that 

is expansion in local co-ordinate at x(P) — ► oo. Its time derivatives (cf. with (|1.19|) ) 

OS k ^ d 2 T 1 

77^ •'• ^ 777-7777777 ( 2 - 2 ) 

j 

form a basis of meromorphic functions with poles at the point x — oo. Suppose now, that these functions are 
globally defined on some curve E and do not have other singularities, except for those at x — oo (the case of a 
single singularity exactly corresponds to the KP hierarchy). If E is sphere (or one considers the dispersionless 
limit of KP hierarchy), this is enough to define all Qk in terms of the " uniformizing" function Qi = A G C. 
Indeed, the set of its powers X k has the same singularities as the set of functions (|2.2() . i.e. these two are related 
by simple linear transformation 

n 1 = x 

d 2 T 



2 = A + 2 



dt\ 



2r o ,2^ (2-3) 



dt\ 2 dhdt 2 



These equalities follow from the comparison of the singular at x = oo part of their expansions in x, following from 
(|2.2|l . Comparing the negative "tails" of the expansion in x of both sides of eq. I|2.3[) expresses the derivatives 
3^ (of S!/c in the l.h.s.) in terms of only those with k — 1 (of A = Sli in the r.h.s.). These relations are 
called the dispersionless KP, or the dKP Hirota equations, and correspond to geometry of rational curve with 
the only marked point (or when the other marked points are simply forgotten), and they can be encoded into 
a "generating form" 

(xi - x 2 ) (l - eDMDMF} = ( D ( Xl ) - D(x 2 )) d tl T (2.4) 

where the operators 

fl W = E (2-5) 

fe>i 

are introduced. 
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In the case of dispersionless Toda hierarchy one should consider sphere with two marked points and corre- 
sponding local co-ordinates z — oo and z — oo correspondingly 1 . Then 



to 
to 

Z — >OQ Z 



dT 1 

dt k z fe +! 



(2.6) 



Considering now two generating functions S = J F zdz and S — J P zdz and, like in (|2.2|l . their time-derivatives 
one gets 



9S h 



Q k = — — = z K 
dt k 



E 

3 

E 



d 2 T 1 

dtkdtj jzi ' 

d 2 ^- 1 

dt k dtj jz3 ' 



fe > 



k> 



(2.7) 



together with 



E 



9 2 J- 1 



»0 = = log Z a + • i 

OTO OtQOtj 1Z 1 

3 J 

= — = log z - 2^ ' 



(2.8) 



dt 



dt dtj jz J 



Suppose again, that <E2|) and are now globally defined meromorphic functions on sphere with two marked 
points. Then it is natural to introduce a uniformizing function w £ C* , using the dipole differential 



dw 



= c?On = —dQr, 



(2.9) 



since the sum of its residues should vanish. It means that 

^3 dtQQtj j z j 



w = ze 

z — >oo 



1 



- E 



3 atndti jz3 



(2.10) 



i.e. w has a simple pole at the point, where z = oo, and a simple zero at z — oo. The Hirota equations for 
dispersionless Toda hierarchy come now from a simple observation that 



<>,, - nun il, r [ - ) I. : 



(2.11) 



for some polynomials P k and P k , similar to (|2.3|) , some of them were explicitly computed in |S] . Formulas (|2.1U|) 
are equivalent to the inverse conformal maps for a simply-connected domain D, which were discussed in the first 
part of this paper pQ. 

For the two-matrix model relations (|1.15|) can be unified into the following expression (a combination of 
formulas (6.21) and (6.3) from pQ) 



log 



w{z) — w(z') 



w(z)w(z') — 1 



2 


1 1 


= log 






z z' 



■ V(z)V(z')^ 



(2.12) 



1 We follow here the conventions of the two-matrix model case 11.21 , corresponding to the independent choice of parameters in 
the vicinity of these two points. 
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which implies an infinite hierarchy of differential equations on the function T . In (|2.12l) we use the operator 

V(zH9 t0 +W f fc\ + ! fc\) (2.13) 

fc>i v / 

discussed in detail in this context in 1 J, whose holomorphic and antiholomorphic parts 

D W = E ^ - D ® = E ^ 9 h , (2.14) 

fc>l k>l 

essentially coincide with (|2.5|) . It is convenient to normalize w(z) by u>(oo) = oo and d z w(oo) to be real, i.e. 

w(z) = -+ 0{\) asz-^oo (2.15) 

where (the real number) r = lim. z _Kx, dz/dw{z) is usually called the (external) conformal radius of the domain 
D or the eigenvalue support of the two- matrix model. Putting z' — > oo in (|2.12|) , one gets 

log \w(z)\ 2 = log \z\ 2 - d ta V(z)T (2.16) 

At z — > oo this equality yields a simple formula for the conformal radius: 

\ogr 2 =d 2 t0 T (2.17) 

Rewriting further (|2.12(l in the form 



(2.18) 



one gets an equality between the holomorphic function of z (in the l.h.s.) and the antiholomorphic (in the r.h.s.) 
function. Therefore, both are equal to a z-independent term which can be found from the limit z — > oo. As a 
result, we obtain the equation 

log ( W [ Z) ; ) = log ( 1 - ~) + D(W (2-19) 



k;(z) - (w(z')) -1 / V z 

which, at z' — > oo, turns into the formula for the conformal map w(z): 

logw(z) = logz - ^9 t 2 ^- a t0 £>(z)^ (2.20) 

where we also used (|2.17ll . Proceeding in a similar way, one can rearrange (|2.19(l in order to write it separately 
for holomorphic and antiholomorphic parts in z': 



- log ( 1 - - L_ ] = D[z)D{z')F 
y w(z)w(z ) I 



(2.21) 
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Writing cqs. (|2.21|l for the pairs of points (z\, z%), {z 2l Z-&) and (2:3, z\) and summing up the exponentials of the 
both sides of each equation one arrives at the relation 

(zi - z 2 )e D ^ D ^ + (z 2 - z 3 )e D ^ D ^ + (z 3 - Zl )e D ^ D ^ = (2.22) 

which is the dispersionless Hirota equation for the dKP hierarchy - a part of the dispersionless two-dimensional 
Toda lattice hierarchy; it is easy to notice, that H2.22|l is just a symmetric form of the equation 12.4JI . This 
equation can be regarded as a degenerate case of the trisecant Fay identity [3] , see formula (|A. 8|) from AppendixlAl 
and discussion of this point in next section. It encodes the algebraic relations between the second order 
derivatives of the function T . As z 3 -> 00, we get these relations in a more explicit but less symmetric 
form: 

x _ e D( Zl )D(z 2 )F = D iZl) ~ Pfo) g T (2.23) 
Z\~Zl 

which makes it clear that the totality of second derivatives Ti$ = dt t dt j T are expressed through the derivatives 
with one of the indices put equal to unity. 

More general equations of the dispersionless Toda hierarchy are obtained in a similar way by combining eqs. 
H2.2L)fl and (|2.21|) include derivatives w.r.t. t and 

{z x - z 2 )e D ^ D ^ = Zl e^ D ( 2 '» f - z 2 e~ d ^ D ^ r 

1 _ e -D{z)D(z)T = J_ e d t0 V(z)F ( 2 - 24 ^ 

ZZ 

These equations allow one to express the second derivatives dt m dt n F, ft m 3f n f with m,n > 1 through the 
derivatives dt dt k T, dtgdj^T. In particular, the dispersionless Toda equation, 

d tl d- tl T = e< r (2.25) 

which follows from the second equation of H2.24(l as z — > 00, expresses dt x dt x T through df a J-. 



2.2 Generalized Hirota equations for the multisupport solution 

To derive equations for the function T in sect. 12. ll we have used the representation 12.21)11 for the conformal 
map w(z) in terms of the derivatives of T and eq. I|2.12l) relating the conformal map to the Green function, 
expressed through the second derivatives of T . In the multiply-connected case the strategy is basically the 
same, with the analog of the conformal map w(z) (or rather of log w(z)) being the embedding of D c into the 
g-dimensional complex torus Jac, the Jacobian of the Schottky double S, see Appendix 1X1 

This embedding is given, up to an overall shift in Jac, by the Abel map IjA.ll) . where the components 
u a (z) — Jg o du> a of the vector (|A.1|I can be thought of as the holomorphic part of the harmonic measure w a , 
discussed in sect. 6 of £Q. The Abel map is represented (see second formula in (6.57) of p^) through the second 
order derivatives of the function T: 



u a {z) — uj a (oo) — / duj a — ~d a D{z)J : (2.26) 
J 00 

and 2 

2Retj Q ,(oo) = m a (co) = —d Ta d a T (2.27) 



2 Following [28|.[T] we are using in this section the time variables t^, which are more strictly defined in general situation, as 
moments adjusted to the proper basis of functions afe(a). In all simple cases with finitely many nontrivial variables they can be 
just identified with the Toda lattice times tj,, for details see the first part of this paper 0. 
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The Green function of the Dirichlet boundary problem can be written in terms of the prime form 
Appendix on the Schottky double: 



6J) (see 

E(z,Q 

Here by £ we mean the (holomorphic) coordinate of the "mirror" point on the Schottky double, i.e. the "mirror" 
of C, under the antiholomorphic involution. The pairs of such mirror points satisfy the condition uj a (£) +uj a (() = 
in the Jacobian (i.e., the sum should be zero modulo the lattice of periods). The prime form 3 is written through 
the Riemann theta functions and the Abel map as follows: 



when the both points are on the upper sheet and 



.Mz)-o,(0) 
h(z) h(C) 



ih(z)h(Cj 



(2.29) 



(2.30) 



when z is on the upper sheet and C is on the lower one (for other cases we define E(z, £) = E(z, Q, E{z, Q = 
E(z, ()). Here 9*(u>) = 65* {u\T) is the Riemann theta function (|A.4|I with the period matrix T a p — 2ni d a dpT 
and any odd characteristics 5*, and 



h\z) = ~z 2 d*A0)dz0J a (z) = z 2 J2 z'k( z )d a d Tk T 

a=l a=l k>l 



(2.31) 



Note that in the l.h.s. of (|2.3U|I the bar means the reflection in the double while in the r.h.s. the bar means 
complex conjugation, the notation is consistent since the local coordinate in the lower sheet is just the complex 
conjugate one. However, one should remember that E(z, Q is not obtained from (|2.29|l by a simple substitution 
of the complex conjugated argument, and on different sheets so defined prime "forms" E are represented by 
different functions. In our normalization (|2.30(l iE(z,z) is real (see also Appendix lA"|l and 



lim 

C^z z 



E(z,C) 



= 1 



(2.32) 



in particular, lirn^oo zE{z, 00) = 1. 

In H2.28fl . the ft,- functions in the prime forms cancel, so the analog of 12.12fl reads 

2 



log 



%{u,{z) - w(0) 



t (w(*) + w(C)) 



log 



V(z)V(C)^ 



(2.33) 



This equation already explains the claim made in the beginning of this section. Indeed, the r.h.s. is the 
generating function for the derivatives Tik while the l.h.s. is expressed through derivatives of the form T a k and 
Tap only, and the expansion in powers of 2, C allows one to express the former through the latter. 

The analogs of eqs. (|zTTd|) . l(TT7|l are 



log 



„(lj{z) - u>(oo)) 



t(u(z) + w(oo)) 



log\z\ 2 + d T0 V(z)T 



(2.34) 



3 Given a Riemann surface with local coordinates 1/z and 1/z we trivialize the bundle of —^-differentials and "redefine" the 
prime form E(z,Q) — > E(z, £ )(dz) 1 / 2 (dQ 1 / 2 so that it becomes a function. However for different coordinate patches (the "upper" 
and "lower" sheets of the Schottky double) one gets different functions, see, for example, formulas 12.291 and 12.301 below. 
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and 

h 2 (oo) 



log 



^*(w(oo)) 

respectively, here vj(z) = 2Reo;(2:) = {w\{z) 1 . . . ,w g (z)) and 



= d 2 TQ T (2.35) 



h 2 (oo) = lim z6* ( [ duj) =-J2 0^ a {Q)d a d Tl T (2.36) 



A simple check shows that the l.h.s. of l|2.35|) can be written as —2 log(z£?(oo, oo)). As is seen from the expansion 
G(z, oo) = — log \z\ — log(iE(oo, do)) + 0(z^ 1 ) as z — ► oo, (iE(oo, oo)) -1 is a natural analog of the conformal 
radius, and (|2.35|) indeed turns to (I2.17|) in the simply-connected case. However, now it provides a nontrivial 
relation on Fap's and !F a iS'- 

(j26*,ad a d Tl Fj [^pe^dpdr^ = ^(^(oo))e<^ (2.37) 

so that the "small phase space" is defined modulo this relation. 

Now we are going to decompose these equalities into holomorphic and antiholomorphic parts, the results are 
conveniently written in terms of the prime form. The counterpart of (|2.19|) is 

log IP^I = log ( 1 - C ~) + D{zMC)T (2.38) 



E((,oo)E(z,0 °V z 
and tending £ — > oo, one gets: 

log S Z '°°? = log z + log E(oo, oo) - d TQ D(z)T (2.39) 
Separating holomorphic and antiholomorphic parts of (|2.38() in £, we get analogs of <|2.21[1 : 

log ^Mb) =log(2 - C)+ ^ (c) ^ (2 - 40) 

-log ^^7'7-j =^MC> (2.41) 

£/(z, O0)£/(00, (J 

Combining these equalities, one is able to obtain the following representations of the prime form itself 



E(z,0 = (z- 1 ~ ^) e -k(D(z)-D(Q)^ 
i£(2,0 = e-^o+^)+S(C)) 2 ^ 



(2.42) 



together with the nice formula 

iE(z,z) ^e-^ 2 ^ (2.43) 

For higher genus Riemann surfaces there are no simple universal relations connecting values of the prime 
forms at different points, which, via l|2.42|l . could be used to generate equations on T . The best available 
relation jS] is the celebrated Fay identity (|A.8(1 . Although it contains not only prime forms but Riemann theta 
functions themselves, it is really a source of closed equations on the quasiclassical tau-function T, since all the 
ingredients are in fact representable in terms of second order derivatives of T in different variables. 
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An analog of the KP version of the Hirota equation H2.22[l for the function T can be obtained by plugging 
eqs. I|2.26[l and 12.4211 into the Fay identity l|A.8() . As a result, one obtains a closed equation which contains 
second order derivatives of the T only (recall that the period matrix in the theta-functions is the matrix of the 
derivatives T a $)- A few equivalent forms of this equation are available: first, shifting Z — > Z — W3 +W4 in 
and putting Z4 = 00, one gets the relation 



(z! - z 2 )e D(zi)D(z2):F 9 / dw+ / dw-Z\ 9 / dw - Z 

\J OO J CO / \J CO / 

+ (z 2 - z^e ^ ^^ 6 (J 2 duj + j 3 duj - Z^j 9 (J ' duj - zj + (2.44) 
+ (z 3 - Zl )e D{z3)D{zi):F 6[ [ * du>+ [ 1 dw - Z ) e( f ' du-Z) =0 



where the vector Z is arbitrary, in particular it can be chosen Z = 0. We see that <|2.22[1 gets "dressed" by the 
theta-factors, and each theta-factor is expressed through for example 

/ dw\ = CX P ( ~^ 2 ^2n a n p d 2 al3 T - 2iti^n a d a D{z)T J (2.45) 

Another form of this equation, obtained from (|2.44|) for a particular choice of Z, reads 

(zi - z 2 )z 3 ~ 1 ei (D(z ^ +D ^ 2 ^h(z 3 )9, f f 1 dw + f 2 du> J + [cyclic per-s of z 1; z 2 , 23] = (2.46) 

\J CO ^ CO / 

Taking the limit 23 — » cxd in 12.44(1 . one gets an analog of (|2.23|) : 



_ fl(C^ + C^-z) gg) -pc^w^ _ 

^dw-Zj^dw-Z) 
D(zA-D{z 2 ) 1 <A 9 9(f Zl duj-Z) 

zi -z 2 zi - z 2 ^ dZ a 0(J u-Z) 



(2.47) 



which also follows from another Fay identity (| A. 9|) . 

Equations on T with ffc-derivatives follow from the general Fay identity (|A.8|) with some points on the 
lower sheet. Besides, many other equations can be derived as various combinations and specializations of the 
ones mentioned above. Altogether, they form an infinite hierarchy of consistent differential equations of a very 
complicated structure which deserves further investigation. The functions T corresponding to different choices 
of independent variables (i.e., to different bases in homology cycles on the Schottky double) provide different 
solutions to this hierarchy 

Let us show how the simplest equation of the hierarchy, the dispersionless Toda equation 12.25|) , is modified 
in the multiply-connected case. Applying d z d^ to both sides of H2.41fl and setting £ = z, we get: 

{dD{z)){dD{z))T = -d z d M log E(z,z) (2.48) 

where dD(z) is the z-derivative of the operator D{z): dD(z) = a' k (z)d Tk . To transform the r.h.s., we use 
the identity (|A.11|) from Appendix 1X1 and specialize it to the particular local parameters on the two sheets: 

\zfd z d s \ogE(z,z) = 6{M %ffl&™^~ Z) + |z| 4 piog0(Z)), a ^a, a (*)%^ (2.49) 
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where 0(Z)) )C(/ g is defined in (|A.12|) . Tending z to oo, we obtain a family of equations (parameterized by an 
arbitrary vector Z) which generalize the dispesrionless Toda equation for the quasiclassical tau-function 



(2.50) 



a, 13=1 



(here we used the z — » oo limits of l|2.2(j[l and H2.43[) h The following two equations correspond to special choices 
of the vector Z: 

a ri 9 ?1 ^+ ]T (lo g 0(O)), a/3 {dv^idpd^F) = e 9 ^ (2.51) 

Q,/3=l 



and 



(2.52) 



a, 0=1 

Finally, let us specify the equation H2.50[l for the genus g = 1 case, in this case there is only one extra variable 
S, which is either Si or IT. The Riemann theta-function 6 (1x7(00) + Z) is then replaced by the Jacobi theta- 
function i9 (dsd To T — Z\T) = $3 (dsd T(l J- — Z\ T), where the elliptic modular parameter is T — 2ni dgJ 7 , and 
the vector Z = Z has only one component. The equation has the form: 



{d S d T0 F + Z\ 2md 2 s T) d z (d s d To T - Z\ 2m d 2 T) ^ 
$ 2 (Z\2md 2 F) 6 T ° 



(2.53) 



-9f logtf 3 (Z\ 2md%F) (d s d T1 F)(d s d fl F) 
Note also that equation l|2.37[) acquires the form 

/ #i(d s d T0 F\2irid 2 s F)\ 2 



(d s d T1 T)(d s d fl T) 



(2.54) 



where = $1 is the only odd Jacobi theta-function. Combining (|2.53(l and 12.54|l one may also write the 
equation 

"tf 3 ( d s d T0 T + Z\ 2m d 2 s T) tf 3 ( dgd T0 J~ - Z\ 2iri d%F) 



d Tl d fl T - 



'tfi (d s d To T\2md 2 T) 
■&[ (0\2md 2 F) 



d 2 (Z\2md 2 s T) 



9| log # 3 (Z\2md 2 s f) 



(2.55) 



d T 

e T ® 



whose form literally reminds (|2.25() . but differs by the nontrivial "coefficient" in the square brackets. In the 
limit T — > ioo the theta-function $ 3 tends to unity, and we obtain the dispersionless Toda equation l|2.25|l . 



2.3 WDVV equations 

Let us finally briefly discuss the Witten-Dijkgraaf-Verlinde-Verlinde or WDVV equations 7 in the context of 
geometry of matrix models. In the most general form [S] they can be written as system of algebraic relations 

TiJ^ x T K = TkT^T!, V I,J,K (2.56) 

for the matrices of third derivatives 

\\FA\jk= mm&TK =FuK (2-57) 
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of quasiclassical tau-function JF(T). Have been appeared first in the context of topological string theories 
[7], they were rediscovered later on in much larger class of physical theories where the exact results can be 
expressed through a single holomorphic function of several complex variables, like we have here for the planar 
matrix models. 

Having the residue formula (|1.16() for the quasiclassical tau-function, the proof of the WDVV equations (|2.56() 
is reduced to solving the system of linear equations [HlEIjj which requires only fulfilling the two conditions: 

• The matching condition 

#(/) = #(«) (2.58) 

and 

• nondegeneracy of the matrix built from l|1.20fl : 

det \\<t>i(x a )\\ ^ (2.59) 



Under these conditions, the structure constants Cf, of the associative algebra 



(2.60) 



responsible for the WDVV equations can be found from the system of linear equations 

<t>l{x a )(j)j{x a ) =22,Ofj(f) K {Xa), Vi„ (2.61) 

K 

with the solution 

C/S = faix^&jixa) {(^Kixa))' 1 (2.62) 

a 

To make it as general, as in jS], one may consider an isomorphic associative algebra (again V x a ) 

(t>l{Xa)(i>j{Xa) = ^ Cfj(C)(f) K {Xa) ' f K) (2.63) 
K 

which instead of l|2.62[l leads to 

CiAO = 2^ Tl — ^ WK{Xa)) (2.64) 

The rest of the proof is consistency of this algebra with relation 

Fijk = Y c ij(0vkl(0 (2.65) 

L 

with 

Vkl{0 = ^ UFkla (2.66) 

A 

expressing structure constants in terms of the third derivatives of free energy and, thus, leading immediately to 
the equations l|2.56|l . It is easy to see that (|2.65|) is satisfied if Tkla are given by residue formula l|1.16f) . which 
can be also represented as 

d z T 1 ^ fdH I dH,,dH K \_ 
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where T a = JYl b ^ a (x a — Xb), see sect. 4.5 of 

Indeed ^UJi requiring only matching #(a) = #(/), one gets 

= E 0j(3; ;^ j / Xq) • (-M^))" 1 • <M^)<M^(* & )r fc (2.68) 

S \ x a) 



K K,a,b 



and finally 



K 



(2.69) 



Hence, for the proof of (|2.56[l one has to adjust the total number of parameters {!/} according to the condition 
l|2.58[) . The number of critical points for the one-matrix model is #(a) = 2n since = in the branching 
points of (|1.8|l . i.e. at = i?(a;) = 0. Thus, one have to take a codimcnsion one subspace in the space of 
all possible parameters of the one-matrix model, a natural choice will be to fix the eldest coefficient of ljl.41) . 
Then the total number of parameters #(/), including the periods S, residue to and the rest of the coefficients 
of the potential will be g + 1 + n = (n — 1) + 1 + n = 2n, i.e. exactly equal to the number of critical points 
#(a). In j^] an explicit check of the WDVV equations has been performed for this choice of parameters, using 
the expansion of free energy computed in |llj . One could try to follow analogously for the two- matrix model 
<|1.2fl with non-hyperelliptic curve (|1.9(l . where the number of critical points in the residue formula #(a) can be 
extracted from pQ, since our analysis here does not use any special properties of the one-matrix case. However, 
the problem is that the number of critical points for the two-matrix model #(a) = 2n 2 + n — 1 (see e.g. formula 
(5.27) from £Q) for generic large n exceeds the naive total number of parameters of the model. 

Another important issue is that formula Ij2.67|l . expressing the third derivatives of quasiclassical tau- function 
T through the quantities 4>i a determined in H1.2()(l . suggests to interpret T as the free energy of a certain 
topological string theory, with three propagators <pi a ending at the same three-vertex: 

'J 



Moreover ^] , one can associate the next term iq of expansion l|1.7f) with the one- loop diagram in this topological 
theory, i.e. with the determinant det/ ja 4>i(x a ): 



Ql,a 

as 

Fx = ~\ log (aV3( x ) ■ det £ - \ log (det (2.70) 

One can therefore conjecture an existence of a diagram technique for calculating the higher genera free energies 
and generating function for the correlators in matrix models. 

3 Matrix models and (p,q)-string theory 

Let us now turn to some physical applications of the geometric picture related to integrability of matrix models, 
which was presented in £Q and developed in the previous section. As a first, already quite nontrivial example, 
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we discuss the correlators in the simplest, so called (p, q) string theory models, directly corresponding to the 
dual matrix model description of two-dimensional gravity |14j . 



3.1 dKP for (p,q)-critical points 

According to existing about fifteen years and very popular hypothesis, see e.g. |15j . the so called (p, q) critical 
points of two-dimensional gravity (corresponding to the generalized minimal conformal matter , interacting 
with two-dimensional Liouville gravity [2]) are most efficiently described by tau- function of the KP hierarchy, 
satisfying string equation (see also sect. 2.3 of £Q). In particular, it means that the correlators on world-sheets 
of spherical topology (the only ones, partially computed by now by means of two-dimensional conformal field 
theory) are governed by quasiclassical tau-function of dKP hierarchy, which is a very particular case of generic 
quasiclassical hierarchy. For each (p, q)-th point one should consider a solution of the p-reduced dKP hierarchy, 
or more strictly, its expansion in the vicinity of nonvanishing t p+q = and vanishing other times, perhaps 
except for t\=x (the so called conformal backgrounds). 

The geometric formulation of result in terms of quasiclassical hierarchy is as follows: 

• For each (p, q)-ih point take a pair of polynomials 

X = \p + ... 

(3.1) 

of degrees p and q respectively. It is convenient to rewrite them as a generating differential 

dS = YdX (3.2) 

• The parameters of quasiclassical hierarchy ("times") and the corresponding "one-point functions" are 
determined by the formulas (|1.14|) . where 

f = X*=A(l + . .. + |?)* (3-3) 

is the local co-ordinate at the point A(Po) = oo. The additional constraint that £ p = £5. or ^_ = (£ p is 
a polynomial in A) means that one deals here with particular p-reduction of the dKP hierarchy. 

• Under the scaling X — > APX, Y — > A q Y, (induced by A — » AA and therefore £ — > A£), the times (|1.14|) 
transform as tk — > A p+q ~ k tk- Then from the second formula of (|1.14|) it follows that the function T scales 
as T — > K 2 ^ p+q ^T, or, for example, as 

Toct'^ f(r k ) (3-4) 
where / is supposed to be a scale-invariant function of corresponding dimcnsionless ratios of the times 

p + g — m 

Tk = tk/t{ +q ~ 1 (|1.14|l . In the simplest (p, q) — (2, 2K — 1) case of dispersionless KdV one also expects a 
natural scaling of the form 

•Foe (t 2K ^) K+ ^ f(t,) (3-5) 
with t; = t2i-i/ (t2K-3)^ K ~' +1 )/ 2 , with the distinguished cosmological time t2K-3 oc A 4 . 

Note also, that the tau-functions of (p, q) and (q,p) theories do not coincide, but are related by the Legendre 
transform ^7], discussed in detail for generic prepotentials in p^. 
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3.2 Pure gravity: (p,q) = (2,3) 

In this case one has only two nontrivial parameters t\ and t 3l and the partition function can be calculated 
explicitly. The times (|1.14|) are here expressed by 

*s = |, h = ^Y 1 -X , t x = jXl-X Y x (3.6) 
in terms of the coefficients of the polynomials (the other coefficients vanish for *2 = ^4 = 0) 

X = X 2 + X 0l F = A 3 + YiA (3.7) 
and, though nonlinear, these quadratic equations can be easily solved for the latter 



1 n~o — — . 1 



X = -^/9t 2 3 -12t 1 -t 3 , Y x = -^tl-l2tx (3.8) 



The second half of residues (|1.14() gives 



— — — A n " — — iiAn 

9*1 8 4 1 
dT 3 



dt 3 -"8 yiX ° 3+ 64 X ° 4 



This results in the following explicit formula for the quasiclassical tau-function 



At t 3 — » 00 (expansion at *i — ► 0) formula (|3.1U|I gives 



(3.9) 



T = 32lo ~ Ut ^' 2 + l 4tl ~ \ Ut " " Tfi (3 - 10) 



T = — ^- 1 + ^ (3.11) 



which is the partition function of the Kontsevich model |18II19) (also identified with the (2, l)-point or topological 

5/2 

gravity). At t\ — > 00 tau function (|3.10(l scales as T cx f/ or partition function of pure two-dimensional gravity 
(14j: expansion at t\ — > 00 gives 



modulo analytic terms. 



3.3 The Yang-Lee model: (p,q) = (2,5) 

The calculation of times according to i|1.14fl gives 



h - ~g^o + -qYsXq - YiXq 

t3 = -jXq - Y 3 X + -Yy 

4 2 3 (3.13) 

t 6 = 7 Y 3 - X 
5 

2 

' 7 =7 
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for the polynomials 

X = A 2 + X a 
Y = A 5 + I3A 3 + YiX 

We see, that <|3.13l) can be easily solved for 



(3.14) 



Y 3 - I (X + i 5 ) 



(3.15) 



and the string equation for Xo is now 

2 

The one-point functions (|1.14|) are given by 



h=-^X*-lt 3 X (3.16) 



dT 1 3 . 3 _ „ 4 9 „„ 3 



— - --YlXo - — x + -r 3 x - --x - -t 3 x 



(3.17) 



where, in the r.h.s.'s, Yj obtained from solving (|3.13|) are substituted. From (|1.15|) . or, differentiating (|3.17|) 
and using H3.16|l . one can easily obtain 4 

d 2 T _ Xo 



dhdt 3 8 u 2 V dt\ 
d 2 F 3 , r3 „ (d 2 T^ 3 



X n J = 3 



dt 3 <9i 3 8 u \dt{ 

However, (|3.16|) cannot be effectively solved explicitly, though it can be solved expanding in t 3 . In the "naive" 
scaling regime (|3.4ll for the partition function of (2, 5) model one gets 

sr^'-i*^ (3 - 19) 

Si 2 ~ 9 1 1 t i/3 J + 9i/ 
Substituting the last expression (using l|3.18|l l into i|3.16l) one gets for the first coefficients of / 

f(r) =fo + fir + \hr 2 + + ^/ 4 r 4 + . . . 

q q o is (3.20) 

f _ y c2/3 f — cl/3 f _ a r2/3 j- _ 10 cl/3 



4 The second equality of 13.181 can be interpreted as dispersionless Hirota equation with an additional KdV (p = 2) reduction 
condition. 
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while /2 (the coefficient at analytic term t\t^ of expansion of T, disappearing from the higher derivatives) 
remains undetermined. The correspondent " scale- invariant" rations are, for example: 

M = Z JL = ^L J^ = G86 fofl = 25 

foh 5' / 2 / 3 125' / 3 / 4 625' flh 14' "' 1 ' ' 

Note, that these scale-invariant rations are given by simple rational numbers, compare to the coefficients (|3.20() 
themselves, which essentially depend upon normalizations of times. In Appendix El it is shown, that expansion 
H3.20|) is in fact related to the generating function in the gravitationally dressed (p, q) = (3, 4) Ising model. 

However, for the non- unitary case one may expect another scaling regime, where the role of times is "ex- 
changed". For the particular case of (2, 5) model this implies that 5 



(3.22) 



^ = 4 /2 f(|,J=4 /2 f(t) 

and string equation l|3.16[l turns into 

t + 5(f") 3 + 3f" = (3.23) 

Solution of (|3.23|) immediately gives for the higher coefficients of expansion f„ = f(")| t of the normalized 
function f the explicit expressions 

f, = ?— 

3 3(1 + 5f|) 

10f 2 

14 : 



9(1 + 5f|)3 (3.24) 
10(-l + 25ff) 
5 27(1 + 5f|) 5 
etc 

through the 2-point function, which satisfies the cubic equation 

3f 2 + 5f| = (3.25) 

while the "normalization" fo and the 1-point function fi remain undetermined from string equation (|3.23|) 
directly. In order to find them, one should use the second equation of (|3.17|l together with the second formula 
of or 

7f - 3tf + 9(f") 5 + 3(f") 3 = (3.26) 

giving rise to 

f _ 9 5 3 3 

f _ 9 f3f 45 4 
fi - - jt 2 f 3 - — f 2 f 3 

In the "phase" with non- vanishing f 2 , this results in 



(3.27) 



, 36 / 3 9 / 3 1 1 / 5 5 r „ 9S x 

fn = \ — , fl = , h — \ — , h = -< h = -\ — ■ f 5 = , •■■ {O.ZOj 

175 V 5' 20' V 5 6' 4V 3' 27' 



5 Such scaling implies also a particular behavior -Fa(*i) = A 7 / 2 z(tiA 3 / 2 ) for the Laplace transformed partition function 
T{tx,ts) = / °° ^e-'3 A F i4 (t 1 ) with "fixed area" (cf. with JM ). 
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Finally, one can easily compute the " invariant ratios" , which do not depend upon normalization of times and 
partition function, for example 6 

fof2 _ _64_ fa fjfs 1 

q ~ 105' f| ~ ' f 2 f 5 ~ s'"- yj 

Analogous formulas for the (2,7) minimal string theory and (3,4) or the Ising model interacting with gravity 
are presented in Appendix iBl 



3.4 Gurevich-Pitaevsky problem and matrix model with logarithmic potential 

In this section we turn to discussion of the Gurevich-Pitaevsky (GP) problem [ST] , appeared originally when 
studying nonlinear waves, and related later to string theory and non-perturbative gauge theories (see e.g. [22]) 
and to the problems of Laplacian growth j2H], see recent paper [23- Solution to the GP problem connects 
different phases dispersionless solutions of quasiclassical hierarchies, corresponding to rational curves £ of 
spherical topology by seweing them via the solutions given by the curves of higher genera. In particular, 
this resolves the singularities of dispersionless solutions. The parameters of solutions appear to be "modulated" 
in accordance with equations of quasiclassical hierarchy. 

Below we consider, first, the original GP problem in KdV hierarchy, whose geometry coincides with that of 
the (2, 5) model of the previous section, and then turn to particular example of two-matrix model l|1.2|l with 
non-polynomial potential. The basic geometric condition for the GP problem is that the finite-gap solutions to 
nonlinear wave equations for the curves of different genera, in order to be correctly sewed when overcoming the 
singularity, must have the same values of the (real parts of the) periods of generating differential § dS over the 
nontrivial cycles. 

In one of two original GP setup's one considers the properties of the cubic string equation Ij3.16|l of the (2, 5) 
model. The corresponding complex manifold arises when one notices, that the polynomials (|3.14|) satisfy the 
equation of (degenerate) (2,5)-curve 

Y 2 = X 5 + (2Y 3 - 5A )A 4 + (Y 2 - 8Y 3 X + 2Y 1 + 10A 2 )A 3 + 
+ (2Y 3 Y 1 - 3Y 2 X - 6YiA - 10A 3 + 12Y 3 X 2 )X 2 + 
+ (Y 2 + 5A 4 - 8r 3 A 3 + 3Y 2 X 2 + mX 2 - 4Y 3 Y 1 X Q )X+ (3.30) 
+2Y 3 Y 1 X% + 2Y 3 X* - A 5 - Y 2 X% - 2Y 1 X* - Y 2 X Q = 
= (A-A )(A-x+) 2 (A-x_) 2 



with 

x ± = -f ± y~\xS-*h (3-31) 

and the sum of the roots, proportional to t$ from (|3.13|) . to be taken vanishing. If all singularities are resolved, 
the curve (|3.30|l would have genus g = 2, and in the case of partial resolution of the singularities it can be 
presented as elliptic curve with an extra " double point" . 

The critical parameters of the solution are determined by vanishing of the periods of generating differential 
dS = YdX. In rational case they can be easily computed, say, using (|3.14() . so that 



A, 

(3.32) 

A; ' V ** / A; 

B Ratios 13,291 coincide with recently computed, using the worldshcet Liouvillc theory, by A. Bclavin and Al. Zamolodchikov, 
see H3. 



jf 3 YdX = yA 7 + ^y 3 A 5 + ^X 3 3 = h 7 + X X 5 + (t 3 + ^X 2 ^j A 
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Figure 1: Vertical section of the matrix model potential (|1.3|) with logarithmic holomorphic part i|3.34[l . which 
has the cusp at z = a (here we took a = 5 and v = 1.) 



where the critical points {A^} are found solving equation Y = 0. Computing l|3.32[l for these limits, one finds 
the critical values _ _ 

h 2 3\/6 y/2 . . 

The ratio of last two values equals to ^^J 27 > found numerically by Gurevich and Pitaevsky [2] , while the first 

one is just the critical value of the equation l|3.16|l . The other two values correspond to x + = x_ and another, 
already nontrivial, case of vanishing of the integral J x x+ YdX — 0. From the point of view of scaling Yang-Lee 
model these critical values correspond to the particular points in parameter space, where the stringy one-loop 
correction to the two-point correlator Xq vanishes. 

Let us now perform a similar procedure for the two-matrix model. Up to now the whole geometric con- 
struction for the two matrix model in was rather implicit, and below we discuss an explicit (and even 
non-polynomial, cf. with discussion in 25 ) example of the potential l|1.3fl with the logarithmic holomorphic 
part 

W{z) = -v (log (l - -) + -) (3.34) 

taking a to be real, see profile at fig. ^ Two its extrema satisfy the equation 

v 

zz — az H — z = (3.35) 

together with the complex conjugated, and are located at z = z = 0, which is the minimum of potential and 
z = z = a — -, which is the saddle point. The corresponding equation of the complex curve (|1.9H for the 
two-matrix model |26j reads 



2 ~2 

Z Z 



(a - -) (z 2 z + zz 2 ) +czz-v (z 2 + z 2 ) + g{z + z) + h = Q{z)z 2 - P{z)z - K{z) = (3.36) 



where last two terms correspond to the points strictly inside the corresponding Newton polygon, see fig. El The 
curve (|3.3()|) is endowed with a standard two-matrix model generating differential 

P 1 




~ zdz =\^ + ^)dz=[— + —^P 2 + 4QK)dz (3.37) 
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Figure 2: Newton polygon (a square!) for the curve (|3.36|) of genus g = 1. 



Rational degeneration of this curve (see general discussion of rational degenerations in Appendix A of PP) can 
be described by conformal map 

u 

z = rw H h v 

r uw (3.38) 
5 = h h v 

w 1 — sw 

where we imply all coefficients to be real; in contrast to the conformal maps considered in Appendix A of 1 , 
the maps (|3.38|l possess singularities not only at w — oo. In the case of rational degeneration of the curve (|3.36|) 
one has 

Vr = A( z - eoV^-e+x^-e-) 

P 1 / P(z) A (3.39) 

where 

A = a+^ (3.40) 
and it corresponds to the following relations between the coefficients in (I3.39|) and l|3.38|l 

£± = V + rs ± 2\fru 



r — u 

to = v + 



(3.41) 



s 

and some more complicated for the coefficients of (|3.36ll . Some useful technical formulas for this case are 
collected in Appendix [U] 

It is easy to see taking the residues of generating differential 13.37f) . that in addition to our common parameter 



tn = </" zdz = res„, zdz H reSoo^-Sdz = res„, zdz res no zdz 

27ri J ]wl=1 2ni qi 2ni + 2ni qi 2ni 92 

,2 



(1 - s 2 ) 2 a 2 + v \ a 2 



(3.42) 



it is convenient to introduce 



11 v 2 

5 = res„. zdz H res Q2 zdz + 2v = a 2 H — 7 — c (3.43) 

2iri 2ni a 2 



Using (|3.43|) and 13.42|1 one can rewrite equation of the curve (|3.36l) as 

- (a - ^ (z 2 z + zz 2 ) + (a 2 + J ^-8^zz-v (z 2 + S 2 ) + 



2 ~2 

Z Z 



(3.44) 
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The last coefficient h is related to the only nontrivial period of (|3.37|) 



S = J- I zdz (3.45) 
2m J a 



In what follows, we fix S = to, corresponding to res q2 zdz = const, then eq. (|3.44|1 acquires the final form 

z 2 z 2 - (a - (z 2 z + zz 2 ) + (a 2 + ^ - t ^j zz - v (z 2 + z 2 ) + t a(z + z) + h = (3.46) 
These formulas give rise to the following system of equations to the parameters of conformal map l|3.38|l 



r + - ^ = As 

1 — s A 



1 - s 2 J V 1 - s 2 



(3.47) 



2 • - . 

r j = t - v 

s z 



/•;/ 



Solving (|3.47|l one gets 



2 V As 

and s 2 = X satisfies the following cubic equation 



2 V As 
l-V As to 



(3.48) 



2A 4 A" 2 2 A 2 A 2 
Eq. (|3.49l) can be solved for t 



^• + ^11-^ + 4 (3.49) 



t = A 2 X 2 ± A 2 X x j [X-\- ( X - 1 + (3.50) 



A / \ A J 

It is easy to see, that two solutions of (|3.5U|) are positive at 



4 +) : < X < 1 y —, X > 1 



A 



r_i 1 2v 2Jv 2Jv 

4 • X<0 > 2~£i~ <X<1 ~A> X>1+ A 



(3.51) 



Thus, the whole picture can be presented as a growing droplet, which is a two-dimensional horizontal surface 
of an "eigenvalue liquid" of the area to filled into the "cup" potential, see fig.H 

It is easy to see, that two phases of real physical interest are: 

• Growing droplet before the "impurance": this corresponds to < X < 1 — and to = ■ The liquid 
of eigenvalues starts to fill the minimum of two-dimensional potential at z — z m j n = and with growth of 
to it level raises till it reaches V^(z sa ddic) = a 2 — tts + 2^ log see fig. |3] In this phase £o > £,+ i.e. (for 
the real choice of all necessary parameters) the cut between £_ and £ + (|3.39|l is located to the left from 
the double point £o- 
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Figure 3: Boundary curves for the first phase. Starting like at the left picture it develops a cusp, when 
approaching z sa ddie = a — ^ (at the pictures a = 5 and v = |.) 
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Figure 5: Numeric values of the integral Re § zdz in different phases: £o > £+ and £o < £- (for a = 5 and 
v = |). It is seen on this picture that they can be easily sewed for different phases, say for Re § zdz = 0. 



• Growing droplet after the " impurance" : this phase corresponds to X > 1 + and to = t\ ' . The liquid 
of eigenvalues fills the impurance at z = a and continues to grow smoothly with growth of to, see fig. 0] 
Now £o < £- i-e. the cut in (|3.39[) is located to the right from the double point. 



*(-) 



Between these two phases it has to pass through the cusp. The cusp point can be "overcome" by the solution 
m spirit of Elj, passing from rational curve out to nonsingular curve (|3.36l) of genus g = 1. The degenerate 
genus g — solutions can be glued through the g = 1 solution only when they have equal values of 



Re <j) zdz 



Re 



Re 



zdz, 



(3.52) 



zdz, £ < £- 



where wir 



are two values of the pullback of the point z = £o to the w-plane. One can see at fig. 



that these values are indeed glued at different phases, in particular when Re / zdz = 0. 



4 Complex geometry of the AdS/CFT correspondence 

Matrix models discussed in U and above present the simplest example of the gauge/string duality |2J. Being 
the " zero-dimensional" example of the gauge theory, they propose a dual description of string theory below two 
dimensions, or the "target-space" theory of two-dimensional gravity on string world-sheets. 

In contrast to matrix model the AdS/CFT conjecture [2Hj deals with the four- dimensional Af — 4 supersym- 
metric Yang-Mills theory (SYM) with the SU (N) gauge group. Again, the main contribution in -^--expansion 
H1.7f) . where closed string loops are suppressed, comes from the planar diagrams when N — > oo at fixed 't Hooft 
coupling X — gyM-N — g s N (coinciding with to for the matrix models, we intensively used in pQ and previous 
sections), while string coupling g s — g YM is equivalent to the quasiclassical parameter H. At X ^> 1 the N = 4 
SYM theory is believed to be dual to string theory in AdS$ x S 5 with the equal, up to a sign (positive for the 
sphere S 5 and negative for the Lobachevsky space AdS^), radii of curvature = A 1 / 4 . Therefore any test of 
the AdS/CFT conjecture implies comparing analytic series at A = (SYM perturbation theory, which is just 
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a direct analog of the counting of diagrams in matrix models as in |14j) with analytic in a' oc worldsheet 
expansion. 

A possible partial way-out from this discrepancy in parameters of expansion can be to consider the classical 
string solutions with large values of integrals of motion (usually referred as "spins" J) on AdSs x S 5 side |3()| . 
whose energies should correspond to anomalous dimensions of "long" operators on gauge side. In this case the 
classical string energy of the form A = yfx£ (^fj^J ma y have an expansion of the form A = J + (7^)' 
over the integer powers of 't Hooft coupling [2], which can be treated as series at X = even at A > 1 
provided large X is suppressed by large value of the integrals of motion 7 J. If it happens (this is not, of course, 
guaranteed) the classical string energy can be tested by direct comparison with perturbative series for the gauge 
theory. 

Such expansion for the classical string energy can be, of course, corrected by quantum corrections on string 
side. Up no now there is no consistent way to quantize string theory in nontrivial background with the Ramond- 
Ramond flux switched on. Therefore, such quasiclassical test of the AdS/CFT correspondence should be treated 
rather carefully. We shall see below, that indeed the quasiclassical picture allows to test the conjecture in the 
first non- vanishing orders, but more detailed analysis requires more knowledge of quantum theory on string 
side. 



4.1 Renormalization and the Bethe anzatz 



The four-dimensional Af = 4 SYM is conformal theory, i.e. /3(<?ym) = 0, but the anomalous dimensions 7 of the 
composite operators, e.g. Tr ($^ . . . ) (x) (where x is some point of the four-dimensional space-time) are still 
renormalized nontrivially. Below we consider the particular scalar operators from this set, though the proposed 
approach can be applied in much more general situation. On string side such operators correspond to the 
string motion in the compact S^-part of ten-dimensional target-space, due to standard Kaluza-Klein argument. 
To simplify the situation maximally, choose two complex $1 = $1 + i$2 and <&2 = $3 + i§4 fields among six 
real 5>j and consider the holomorphic operators 



Tr ($i$i$i$ 2 *2*l*2*2*l*l*l*2 ) (x) 



(4.1) 



which can be conveniently labeled by arrows as ITTT4ITI4TTT4 • ■ ■) G (C 2 ) 

These operators do not have singularities at coinciding arguments of all fields, since the kinetic terms 
for the holomorphic fields do always have the $$-structure. The holomorphic subsector is "closed" under 
renormalization and the anomalous dimensions are eigenvalues of the 2 L x 2 L mixing matrix 



H = 



X 



16tt 2 



0(A 2 ) 



(4.2) 



1=1 



which is, up to addition of a constant, the permutation operator in (C 2 ) <g> (C 2 ) , whose appearance is determined 
by structure of the <f> 4 -vertex in SYM Lagrangian, or the Hamiltonian for Heisenberg magnetic 35 . 

It is well-known, that the matrix (|4.2|l can be diagonalized using the Bethe anzatz 0B], (see e.g. [37] for 
present status of this technique and comprehensive list of references) . Integrable structure of Heisenberg spin 
chain is encoded in the transfer matrix 



f(u) = Tr 



C7l 



(4.3) 



'On integrability on the string side of duality see, e.g. 1321 . 
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which can be presented as an operator-valued matrix in the auxiliary two-dimensional space 

/ A(u) B(u) \ 

T(u) =TV =A(u) + D(u) (4.4) 

\ C(u) D(u) j 

with four "generators of Yangian" A(u), B(u), C(u) and D(u), satisfying the quadratic RTT-algebra. Expansion 
of the operator l|4.3|l in spectral parameter at u = | generates the local charges 



L 

i ' 



T(u) = ( u + - ) Uexp 



n—1 ^ 



(4.5) 



where U = e lP is the shift operator (by one site along the chain) and the Hamiltonian (|4.2|) Qi oc H is one of 
the local charges. Contrarily, expansion of (|4.3(l at u — oo produces the non-local Yangian charges, whose role 
in the context of N = 4 S YM was discussed in |3Sj ■ 

Bethe ansatz diagonalizes the whole spectral-parameter dependent operator (|4.3|l . I|4.5|l 

f(u)\u 1 ...uj)=T(u)\u 1 ...u J ), (4.6) 
explicitly constructing the eigenvectors \ui . . . uj) in terms of the operators l|4. 4f> 

\ u i ■ ■ ■ u .j) = B( u i) ■ ■ ■ B(uj) ITTTTT • • •) (4.7) 

acting on the ferromagnetic vacuum ITTTTT ■••)• Any Bethe vector l|4.7|l depends on J rapidities, and the 
corresponding eigenvalues of the operator (|4.6|l are 

Now since any eigenvalue T(u) of the operator l|4.3|l must be a polynomial in u, the cancelation of poles at 
u = Uj in the expression l|4.8|l gives J conditions: 

u j + i/2 \ L = " ur-Uk + i 

which are the Bethe equations. 

Hence, the eigenvectors of Hamiltonian H4.2J) are parameterized by Bethe roots {u\, . . . ,uj}, (J < \L due to 
an obvious Z^-symmetry, in particular permuting two ferromagnetic vacua |T I I l'T ■ ■ ■) ||||44 • • •)), satisfying 
H4.9fl . After taking the logarithm, eqs. I|4.9|l acquire the form 

fu.j+i/2\ „ , Uj—Uk + i , 

L log ^ L- = 2m nj + V log -2 k - : 4.10 

y Uj -i/2j ^ uj-u k -i 

with appeared mode numbers nj G Z being some integers, never vanishing simultaneously for a nontrivial 
solution. For the diagonalization of the mixing matrix of the operators (|4.1|) in SYM theory, equations (|4.1U|) 
are supplied by the "trace condition" 



nu-j + i/2 
j=i 



i.e. integrality of the total momentum of magnons: ^-eZ 
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Figure 6: Riemann surface S, which is a double-cover of the s-plane cut along the Bethe strings (four slightly 
curved lines on each sheet in this example), which cross the real axis at x = o^T- There can be also possible 
cuts along the imaginary axis, which can be however eliminated after appropriate choice of the homology basis 
in see [33] . The upper sheet is physical for the Abelian integral G(x), while on the lower sheet the 

resolvent has extra singularities. 

The energy of the eigenstate of the Hamiltonian 1)4. 7)1 is given by sum over magnons, i.e. oc ^2j=i dt 

for the quasimomentum expip(w) = ' anc ^ ^ means that the anomalous dimensions for the mixing matrix 
14.2)1 in the first order of the SYM perturbation theory can be written as 

i = &t<rTyi + oW (4.12) 

j=l J 

and these 7's one needs for testing the quasiclassical AdS/CFT conjecture. 

For comparison with dual (classical) string theory we are interested in the long operators with L — > 00, for 
which it is known empirically that the Bethe roots are typically of the order of uj ~ L. Rescaling Uj = Lxj, 
and omitting higher in j- terms, one gets from 1 )4.101 the equations of more simple form 

1=2^ + 1^^4- (4.13) 

J kj=J J 

Like in the matrix model case (see detailed discussion in p^), if the second "interaction term" in the r.h.s. is 
absent, the solution to 1)4.13)1 is given by Xj — for each rij (filling the extremal), and when one switches 
on the interaction, the roots corresponding to rij (if there are several for each value) will " concentrate" around 
2"^— forming the so called " Bethe strings" , shown at fig. HJ In other words, if several different Xj 1 Xj m have 
the same mode number rij = = . . . = rij m , these m Bethe roots combine into a Bethe string with the mode 
number rij m 0. 

Introduce at L — > 00, as in sect. 3 of 0, the density 

1 ^ f J 



i=i 



2G 



or resolvent 



= I £ - J_ = / , -L / ^ = ' (4.15) 

Lj^x — Xj J c x - £ 2m J c L 



3 = 1 



and consider the finite number of different Bethe strings nj 7^ n;/, with = 1, . . . ,K despite the infinitely 
many roots J. Then in the scaling limit (J — > 00 with fixed finite K) the total eigenvalue support is C = 
Ci |J . . . 1J Ck , where on each component one gets from H4.13(l 

2 f *t£M = G( x+) + G(x-) = --2mi u 16C, (4.16) 
J c x - £ x 

where G(x±) are values of the resolvent on two different sides of the cut along the Bethe string. 

The integrality of total momentum condition, using the Bethe equations (|4.13|) . acquires the form 



1 ^ 1 . \ r 

-V" — = 27ry^ni / p(x)dx = 2mn, meZ (4.17) 
1 / G(x)dx 



3=1 3 1=1 

or 

2wm, ni,m£Z (4-18) 



2ni 



c 



Different ni ny on different parts of support C; f] Ci> = mean that, in contrast to the matrix model case, 
discussed in detail in sect. 3 of pQ, G(x) = J dG is not already a single- valued function, but an Abelian integral 
on some hyperelliptic curve E, introduced by hands, according to the number of nontrivial Bethe strings 

2K 

y 2 = R 2 k(x) = x 2K + rn 2 ^ 1 +...+r 2 K = Y[{x- xj) (4.19) 

3=1 

It means, that resolvent G(x) in contrast to the case of one-matrix model does not satisfy here a sensible 
algebraic equation, but can be still determined from its geometric properties. 

Equations (|4.15|) . (|4.1ti|) and (|4.18|l can be solved after reformulating them as a set of properties of the 
meromorphic differential dG of the resolvent (|4.15|) : 

• The differential dG is the second-kind Abelian differential with the only second-order pole at the point 
Po, (x(Pq) — on unphysical sheet of the Riemann surface E, see fig. EJ), in particular it means that 
§ A dG = for all ^4-cycles, surrounding the Bethe strings. The absence of poles on physical sheets is 
determined by properties of the resolvent, but on unphysical sheets then any differential with vanishing 
^-periods must have singularities; 

• In addition to vanishing A-periods, the differential dG has integral _B-periods 

dG = 27r(n 4 - n K ) (4.20) 

'Bi 

More exactly one can write [3"3"| 



dG = 2m j , j = l,...,K + l (4.21) 

where B'j is the contour from oo + on the upper sheet to oo_ on the lower sheet, passing through the j-th 
cut, so that Bj — B'j — B' K , for j = 1, . . . , K , see fig. 
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dG has the following behavior at infinity 



dG = — + ( 4 - 22 ) 

x^oo L X 

and the Abelian integral G(x) itself is fixed by (|4.18|l 

G{x) = 2nm + [ dG, or G(0) = 27rm (4.23) 
Jo 

The general solution [33] for the differential dG on hyperelliptic curve l|4.19[) . satisfying the above requirements 
H4.20f> . I|4.21|l . (|4.22() and (|4.23() is read from a standard formula for the second-kind Abelian differential on 
hyperelliptic curve 

K—l j. i 



dG = -^ s (l- ^ ) + r -^=L- + y a k X —^ (4.24) 

together with the extra conditions, ensuring, in particular, the single- valuedness of the resolvent on "upper" 
physical sheet 

dG = 0, i = 1,2,..., K- 1 (4.25) 

Ai 

which is a system of linear equations, to be easily solved for the coefficients {afc}. The rest of parameters is 
" eaten by" fractions of roots on particular pieces of support 



[ p(x)dx = (f xdG, 



l Cj 2tti J Aj (4.26) 

j = l,...,g = K -I 



the total amount of Bethe roots (|4.14ll . and the total momentum (|4.18|l . 

The energy or one-loop anomalous dimension for generic finite-gap solution |33j can be read from H4.12|l . 

X I dx nt-.\ X ( r 2K-2 r 2K-l a-i 



Btt 2 L J c 2nix 2 8Tr 2 L\4r 2K 16r| K ) 

The anomalous dimensions defined by 14.27|) are functions of the coefficients of the embedding equation (|4.19|) 
and ai which again is expressed through these coefficients by means of l|4.25(l . The moduli of the curve (|4.19() 
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are themselves (implicitly) expressed through the mode numbers rij and root fractions Sj via l|4.2U[l or (|4.21l) 
and (|4.26() (together with the total momentum (|4.18l) and the total number of Bethe roots l|4.15|l ). 

It is interesting to point our, that formula l|4.24ll is a particular case of the two-point resolvent for the planar 
matrix model Wo(x,£ l ), or so called Bergman kernel u(P, P') = d P d' P log E(P, P') (see also formula <|A~T3)> from 
Appendix 1X1). which on a (tensor square of) the hyperelliptic curve 14.19f) reads 



dxd^ 



2(x-0 2 2y(0 
dxd£ 



1 



1 



In 



a = l 



W Q (x,£)dxd£ = 

1 

(f - x)(x - X Q ) 



R(x) 



n—1 „ 



dx 



\ dxd£ R'{x) 



' Ai 
2n n-1 



(x - x Q ) 2 y(a;) 



dxd£ 



(4.28) 



-y(x)dx^2^2dwi(£) f 



dx 



ct—1 i—1 



a, - y- a ) 2 y{x) 



The detailed discussion of this issue can be found in ^3] ■ 



4.2 Geometry of integrable classical strings 

Formulas H4.24JI - i|4.26J) of the previous paragraph show, that the general solution for anomalous dimension 
of long operators l|4.27|l is expressed through the integrals of motion on some classical configurations of the 
Heisenberg magnet 8 . In the dual string picture one has the classical trajectories of string, moving in (some 
subspace of) AdS§ x S 5 and the finite gap solutions to string sigma-model in the Lobachevsky-like spaces were 
first constructed in |4(J|. Being slightly modified, it can be easily applied to the case of compact S 20-1 sigma- 
models. In this section, following we are going to show how these classical solutions can be compared 

with the quasiclassical solutions on the gauge side. 

In particular subsector of only two holomorphic fields one gets the S 3 C S 5 sigma-model (in the AdS^-sector 
the only nontrivial string co-ordinate on the solution is "time" Xq = t) with the action 

S =1^J dadT [(^) 2 - (daXof] , (4-29) 
{^2,Xf — 1). Provided by identifications (since S 3 is the group-manifold of the Lie group SU(2)) 

e 517(2), 



X x + iX 2 X 3 + iX 4 \ Zi z 2 

-X 3 + iX 4 X x -iX 2 j ~ { -z 2 z x 



1 / Z x dZ x + Z 2 dZ 2 Z x dZ 2 - Z 2 dZ x \ 
J = g dg = - - - <= su(2). 

\ Z 2 dZ x — Z x dZ 2 Z 2 dZ 2 + Z x dZ x j 

it is equivalent to the SU (2) principal chiral field with the Lax pair, 

A iS± • a 
J+ix) = — = 

d+J- -d-J+ + [J+,J-] = 
d+J- +8-J+ = 



(4.30) 



(4.31) 



which has two simple poles at values of string spectral parameter X — X(P±) — In different words 

8 Such correspondence with classical solutions was first noticed in 1391 for the non-linear Schrodinger equation. 



29 




Figure 8: Riemann surface T, which is a double cover of E with a single cut. It has an involution exchanging 
S <-> E. 

such sigma-model is equivalent to a system of iwo interacting relativistic spins S + and S_: 

2A 

<9+S_ + -pS_xS + =0, 

2A ( 4 ' 32 ) 
<9_S+ - — S_ x S+ = 0. 

which in the "non-relativistic limit" 

S±=S 4 SXS -^^ S + - (4 - 33) 

degenerates into the Heisenberg magnet, similar to the procedure, studied in the papers |42|. 

However, this method can be used only for the group-manifolds. Nevertheless, in general situation the string 
sigma-model solution for the complex co-ordinates Zj(t, a) and Zj(t, a) on S 20 ^ 1 , (constraint by J^i \Zi\ 2 = 1) 



Z I (a ± )=ri?(q I ,a±) Z x (tr ± ) = r/Tfe a±) = nT( 9j , a±), 7=1,...,£> (4.34) 
(<r± = i(cr ± r)) can be found 00] in terms of the Baker- Akhiezer (BA) functions 

T(P,a ± ) p = p e fe±ff ± + cxe°+( p ^++ n -( p ) CT -0(u;(P)+V +( j + +V_ ( 7_) (4.35) 

where, w(P) is the Abel map (|A.1(1 and 6* (<*>) is the Riemann theta function l|A.2|) (see Appendix lAl for details), 
defined on double cover V (branched at P+ and P_ ) of a Riemann surface E (see fig. |HJ . For only two complex 
co-ordinates Zi (like in the S 3 case) the curve E is hyperelliptic and directly related with the curve l|4.19|l of 
the Heisenberg chain. 

More strictly, the string hyperelliptic curve E turns into the spin chain hyperelliptic curve (|4.19|l . depicted 
at fig.Elonly in the limit A/A 2 — > 0, consistently with (|4.33() . Schematically it can be drawn as on fig. 03 which 
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Figure 9: Riemann surface E for classical string theory. Compare to fig. the pole on unphysical sheet is 
replaced by a cut between P + and P_ of the length ^5 . Two copies of this Riemann surface, glued along this 
extra cut give rise to the full string curve T from fig. [S] 



is quite similar to the spin chain curve from ng.[S]except for the pole on unphysical sheet, which is now replaced 
by an extra cut of the length ^-J- 

The BA function H4.35[l (and hence the solution to sigma-model) is constructed in terms of two second-kind 
Abelian differentials dQ± on Riemann surface T 

dfl± = ±dk ± (1 + 0(fc± 2 )) , / dfl± = (4.36) 

P^ p ± J A 

with the only second-order pole at P± respectively; V± = ^ B dil± are the vectors of their P-periods. 

The proof of the fact that formulas H4.34J) are solutions to the sigma-model, satisfying classical Virasoro 
constraints, is based on existence of the third-kind Abelian differential dfl on E with the simple poles at P± 
and zeroes in the poles of the BA functions T and conjugated T. Then one may define the string resolvent or 
quasimomentum by the following formula 

dQ = - (dn + -dn_) 

1 r (4-37) 

dQ = -dk± (1 + 0(fc± 2 )) , i dQ = 
P—>P± 2 J A 

For periodic in a solution, as follows from (|4.35() , the P-periods of the resolvent ^ §-^dQ gZ are integer- valued, 
and for the periodic solutions one can write 

dQ = 2Z-dQ (4.38) 

where brackets mean the average over the period in cr-variable. 

The BA function (|4.35|) satisfies the second-order differential equation 

{d+d- +U)?(P,(t±) = 0, Per (4.39) 
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where U oc J^i (d+Zid-Zi + d-Zid+Zi) . This fact and the Virasoro constraints X J^i \d±Zi\ 2 = A are guar- 
anteed by the properties of the differential l|4.38() and existence of the function E on £ with D simple poles 
{qi} and the following behavior at the vicinities of the points P±: E = E± ± + . . ., (cf. with |4(J) ) . 

The normalization factors in the expressions for the sigma-model co-ordinates <|4.34[1 are determined by the 
formulas 

9 res„. EdVt „ 
r i= E __ E » I=l,-.,D (4.40) 

where E± = E(P±) and normalizations (|4.40() satisfy X^iLi r] = 1 due to vanishing of the total sum over the 
residues res (EdQ) = 0. 

Rescaling d± ^ d T ± d a , T ^ e^ T *, T -> e~^ T ^, U -> J7 - in the limit A/A 2 < 1, one gets 
from H4.39f) the non-stationary Schrodinger equation 



(d T - dl + U) * = 

(-flr - a 2 + c/) * = o 



where now both BA functions W and ^ can be defined on Riemann surface £ (see fig.E|), with the ends of the 
extra cut P± are shrinked to a single point Pq, with the expansion at the vicinity of this point (with new local 
parameter fc(Po) — oo) 



p^p \ k k 2 



p->p \ k k 2 

Substituting expansions 14.42|) into (14.41ft one gets 

Oa Oa 



(4.42) 



^-2^-1^ + ^ = (4-43) 
or Oa da 



2 

— 7T~ + 2 15 S~T" + = 

or oa Oa 2 

For the conjugated functions one has tp\ + ipi = and then ^ (-02 + V"2 + V'iV'i) = — Therefore, one gets 
an expansion 

= 1 + -^ + ... (4.44) 

P^Po fc 2 

The differential l|4.38|l in this limit turns into dQ = dG and function E acquires a simple pole at Po, i.e. 
E 1 = k + . . ., if written in terms of new local parameter fc. From vanishing of the sum of the residues of 

P— »Pn 



differential Edfl one gets now 



U = res Po Edfl = - ^ E—^——dG cx ^ (4.45) 



(**> 

91 x ' I 



It means, that ^j(r, cr) cx ^(g/) satisfy some vector non-linear Schrodinger equation 03] 

(dr-dl + Y^^J? j */ = (4.46) 
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In the case of D = 2 the curve £ is hyperelliptic and one can take the function E with the only two poles, see 
below. Then l|4.46[) turns into the ordinary non-linear Schrodinger equation, which can be transformed to the 
Heisenberg magnetic chain [441 14*5] : 

l*| 2 °c $1 

- (4.47) 
- I'd* oc (S CT • S x S CTCT ) v ' 

and so on, which is a gauge transformation for the Lax operators. 

An equivalent way to describe the classical string geometry was proposed in |33| and was based on reformulat- 
ing of geometric data of the principal chiral field l|4.31|> in terms of some Riemann-Hilbert problem. The spectral 
problem on string side (a direct analog of the formulas (|4.15|l . (|4.16|l . (|4.18() and (|4.27(l 'l can be formulated in 
the following way. Let X and Q(X) be string spectral parameter and resolvent, equal to the quasimomentum of 
the classical solution (maybe up to an exact one- form). The spectral Riemann-Hilbert problem on string side 
can be written as [23 

— i g{x)dx = i. + ^-^ 

2m Jc A 2A 

1 £ dXG(X) 



and 



JW/ C X 2 ™ ^ 

2tdx g{x) _ 

g(X + ) + g(X^)- 2™, = -£— t (4.49) 

where we introduced the notation t = . c 4 A ■> . 

Comparison of the string (|4.48|l . (|4.49|) and gauge l|4.16[l . Ij4.18(l Riemann-Hilbert problems implies that 
x = X + y together with G(x) = Q{X) is an exact change of variables for the quasiclassical AdS/CFT 
correspondence. This in particular automatically implies the involution X <-> -L on the string theory curves T, 
see fig. El Indeed, in terms of k± one has 

X = ±Vi + — + ... (4.50) 

p^p± k± y ' 

then 

x = X + ± p = p± ±2V-t±± t] L + ... (4.51) 

i.e. the function E = — (up to an overall constant) satisfies all desired properties for the modified construction 
from jlU] we presented above, e.g. when the cut between P + and P_ on fig. [5] shrinks to a point Pq on fig. El 
with x(Pq) = 0, the function E acquires a simple pole at this point. 

Proceeding further one gets 

dX 





dX 












~ ~x 





= Vx 2 - At (4.52) 



and combination of the first and third lines in l|4.48|l gives 



— / g(X)dX = - + ^± = - + tt (4.53) 

2mJ K ' A 2A A / 2niX 2 K ' 

or 

± I dxG(x) = L (4.54) 
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The second line of 1(4.48(1 is then 

1 / dxG(x) . . 

* = 2nm (4.55) 



2™ J yjx 2 - At 

where the integral is taken around the cut between the points —2\fi and 2\fi in the x-plane, and the third line 
of 14.48( 1 gives 

I dxG ^f x l)=A-L (4.56) 



J 2m \y/x 2 -At 
The "string Bethe" equation on the cuts (|4.49|) turns now into 



G{x + ) + G{x_)-2imi = -= : (4.57) 

\/x 2 - At 

We now see from l|4.54(l . (|4.55(l . ((4.56|) and (|4.57(l that the classical string theory spectral problem indeed is 
identical to the quasiclassical Bethe equations on gauge side upon replacements 

1 1 _ 1 2t 

x ^x 1 - At x x 3 

L^A ( 4 - 58 ) 
7^ A-L 

In other words, this leads to a nonlinear relation 

A — L = A) (4.59) 

where T(X,L) — 7 + 0(X 2 ) should be compared with the first perturbative contributions to the anomalous 
dimension of the supersymmetric gauge theory. 

A simplest non-trivial example of such relation is the solitonic limit of small number of Bethe roots, leading 
to the "modified" BMN formula [33] 



A-L = E^(V 1 + ^- 1 ) (4 - 60) 

for J = ^ fc Nk expressed as a total amount of "positive" nk > and "negative" rik < massive oscillators |46| . 
Formulas (|4.59() and ((4.60(1 show, that the solution for A of classical string theory is given in terms of the highly 
non-linear formulas, and the oscillator language of (461 IT7) is rather an effective tool for description of certain 
quasiclassical modes of an integrable string model in pp-wave geometry, than an exact world-sheet quantization 
of the theory. 



4.3 Beyond the SU(2) subsector: overparameterized curves 

It turns out to be very hard to write down an algebraic equation for the curve beyond the SU (2) subsector. 
Up to now, the only achievement in this direction was related with the highly "overparameterized" curves, 
i.e. determined by equations, giving rise to much higher genera and much larger set of parameters, which is 
necessary for solving the problem. These curves should be further constrained in order to fit with the moduli 
space of necessary dimension for the desired solution. The situation here is a bit similar to the two-matrix 
model, see discussion of the complex curve of two- matrix model in . 

The simplest example of such overdetermined situation is given by the so called Baxter curves, which are 
determined by the Baxter equation 

e lp + e- lp = Tr Q(x) (4.61) 
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Figure 10: Newton polygon for the non-degenerate curve, defined by the equation (|4.63|1 . The smooth genus g 
given by the number of integer points inside the polygon is g = 3(2K — 1) = 6K — 3. 



for the quasimomentum p = G — w-, which define the finite genus curve only for the chains with finite number 
of sites L, this was used for similar purposes in 03- Otherwise, if L — > oo, l|4.61|l corresponds literally to a 
complex curves of an infinite genus, and cannot be used effectively for any of practical purposes. Differently, 
the algebraic function p has essential singularities on the "minimal" curves, like (|4.19|) . and it is not generally 
algebraic, unless the corresponding equation gives rise to an overparameterized curve of an infinite genus, where 
the essential singularity if " dissolved" . 

Moreover, instead of 14.19JI one can write down an equation satisfied by the meromorphic functions p' = x 2 ^ , 
or 4^ for the resolvent l|4.24|l . Clearly, it is has the form 

P 2 (x)(p') 2 + P (x) = (4.62) 

where P%{x) = B,2k{x) and Po(x) are polynomials of degree 2K . Note, however, that equation l|4.62|l corresponds 
to a curve of genus g = 2K — 1, instead of g = K — 1 for the minimal curve 14.19fl . i.e. it has K extra 
"false" handles with corresponding extra parameters. We shall refer to the curves of the type (|4.ti2|) as to 
"overparameterized" curves, being "in between" the minimal and the Baxter ones. 

Beyond the SU(2) case the overdetermined curves of the type l|4.62|l were proposed in |3SJ m the form of 
quartic equation 9 

Pi{x)(p') A + P2(x)(p') 2 + Pi(x)p' + P (x) = (4.63) 

with all polynomials Pi(x), i = 1, . . . , 4 of degree 2K. It is easy to see from the corresponding Newton polygon, 
see fig-EH that the genus of H4.63[) is g = 6A" — 3. 

If one additionally requires, however, the only 2K branch points for the curve, defined by quartic equation 
H4.63fl . the Riemann-Hurwitz formula 2 — 2g = 1(2 — 2go) — j^B.P. gives 

2 - 2o = 4(2 - 0) - 2K = 8-2K, or 
9 = K - 3 

i.e. one has to constraint g — g = 5K parameters. In 49: these constrains are chosen in the form, suggested by 
comparison of (|4.62(l with the minimal curve Ij4.19|l . 

9 A similar analysis for the most general case of octic equation can be found in lolll . 
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5 Conclusion 



In this paper, as in its first part , we have concentrated mostly on the geometric properties of the quasiclassical 
matrix models, closely related with their integrability. The corresponding geometric data are encoded in one- 
dimensional complex manifold - a complex curve or Riemann surface, endowed with a meromorphic generating 
one-form. The free energy of matrix model is given by a prepotential, or quasiclassical tau-function, completely 
determined by these geometric data. The geometric definition of the quasiclassical tau-function satisfies the 
simplest possible integrability relations - the symmetricity of its second derivatives, guaranteed by the Riemann 
bilinear relations for the meromorphic differentials on this curve. The derivatives of prepotential can be express 
in terms of residues, periods and other invariant analytic structures, and as a consequence, the quasiclassical 
free energy of matrix models satisfies certain non-linear integrable differential equations, which were studied in 
the paper. 

In this respect an immediate question arises, concerning the rest of the expansion 1)1.7(1 : is there any natural 
geometric picture for the correlators of matrix models beyond the planar limit? It does not have any clear 
answer at the moment, except maybe for the string one-loop correction F%, is given by certain determinant 
formulas, as l|2.70[l . which in certain sense describe the fluctuations of the quasiclassical geometry, see also 
[511 I13| . As for the whole sum l|1.7|) , it is rather described by a quantization of the presented above geometric 
picture, and some step towards this direction were discussed, say, in |52) . From the point of view of integrable 
systems quantization of geometry of one-matrix model corresponds to the so called stationary problem in the 
Toda chain hierarchy |53| . while for the two- matrix model case one has to consider the same problem in the 
context of the full hierarchy for the two-dimensional Toda lattice. 

Even more questions arise concerning the quantum generalization of quasiclassical picture of the AdS / CFT 
correspondence (see discussion of some of these issues in 54 ). While one can certainly move in this direction 
similarly to studying matrix models within the 1/iV-expansion (here the same role in some sense is played 
by the finite-size corrections to the quasiclassical solutions of Bethe equations), the understanding of the full 
quantum picture necessarily requires at least some progress in formulation of quantum string theory in non- 
trivial backgrounds. We believe that such progress is possible, but this is certainly beyond the scope of particular 
questions, discussed in this paper. 
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discussed in this paper. This work was partially supported by RFBR grant 04-01-00642, the grant for support 
of Scientific Schools 1578.2003.2, the NWO project 047.017.015, the ANR-05-BLAN-0029-01 project "Geometry 
and Integrability in Mathematical Physics" and the Russian Science Support Foundation. 

Appendix 

A Theta functions and Fay identities 

In this Appendix we present some definitions and useful formulas for the Riemann theta-functions, basically 
taken from the book 5 . One has first to consider the embedding of the Riemann surface E into the g-dimensional 
complex torus Jac, the Jacobi variety of E. This embedding is given, up to an overall shift in Jac, by the Abel 
map P i— ► <j}{P) — (uji(P), . . . ,uj g (P)) where 




Po 



(A.I) 
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The Riemann theta function 9(u>) = 6(u>\T) is defined by the Fourier series 

V\f>) = }^ e (A.2) 
with the positively definite imaginary part of the period matrix of S 

lvaT a p oc / du) a A dujp > (A. 3) 

The theta function with a (half-integer) characteristics 5 = (Si, 82), where 5 a = T a p5i t + S2, a and Si, 82 € ^Z 9 
reads 

9g(0j) = ^di-T-Si+^i-^+Si)^^ + (5) = ^ e ^(n+S 1 )-T-(n+S 1 )+27ri(n+S 1 )-(UJ+S 2 ) ^ 

n£Zs 

Under shifts by a period of the lattice, it transforms according to 

5 (u, + e a )=e 2 ^0 5 (u,) 
+ lapep) = e og(u;) 

The prime form E(z, £) is defined as 

P ( P > P ) = . : = /_ : . = (A.6) 



where 0* is any odd theta function, i.e., the theta function with any odd characteristic S* (the characteristics is 
odd if 481 • #2 = odd). The prime form does not depend on the particular choice of the odd characteristics. In 
the denominator of 1A.6I) we used the notation 



. -KM = "- M 



duj a 

for the set of ^-constants. 



(A.7) 



The data we use in the main text contain also a distinguished coordinates on a Riemann surface: the 
holomorphic co-ordinates z and z on two different sheets of the Schottky double (see details in pQ), and we do not 
distinguish, unless it is necessary between the prime form (|A~6)l and & function E(z,() = E(z,C,)(dz) 1 / 2 (dC,) 1/2 
"normalized" onto the differentials of distinguished co-ordinate. 

Let us now list the Fay identities [5] used in |2Hj and above in the paper. The basic one is the trisecant 
identity (equation (45) from p. 34 of 0) 

9(lji - u> 3 - Z) 0(u> 2 - w 4 - Z) E(zi, z 4 )E(z 3 , z 2 ) 
+ 6(u>i - u> 4 - Z) 0(u> 2 - u> 3 - Z) E(z u z 3 )E(z 2 , z 4 ) (A.8) 
= 9(oji + uj 2 - ^3 - w 4 - Z) 0(Z) £(zi, z 2 )E(z 3 , z A ) 

where Wj = 07(2^). This identity holds for any four points Zi, . . . , Z4, on a complex curve and any vector Z G Jac 
in Jacobian. In the limit z 3 — > 2:4 = 00 one gets (formula (38) from p. 25 of [S]) 

0tCdw+ f* 2 du,-Z)0(Z) i?^,^) 



0(£ du> - Z) 0(j^ 2 du - Z) £7(2 1|0 o)£(«2,oo) 
drt^(oc) + £ <M«>) d Za log J ~ — 



(A.9) 



a=l 



<?(Cdw-Z) 
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where 

dn&>*>\oo) = 4 log (A. 10) 

E(Z,Z2) 

is the normalized Abelian differential of the third kind with simple poles at the points with co-ordinates z\ and 
Z2 and residues ±1. 

Another relation from we used above (see e.g. (29) on p. 20 and (39) on p. 26) is 

P(Z)E 2 (z l Z2 ) -= w (*i.*a)+ E flog0(Z)) iO p<M*i)<M*O (A.ll) 

where 

<*"<*».■* 

and 

u{z 1 ,z 2 ) = d Zl d Z2 log£ , (zi,z 2 ) (A. 13) 

is the canonical bi-differential of the second kind with the double pole at z\ = z 2 or the Bergman kernel. 

B Correlation functions in (2,7) and (3,4) minimal string models 

(2, 7) model. The (p, g) = (2, 7) model is not much different from the Yang-Lee case of (2, 5) theory considered 
in sect. 13.31 The polynomials (|3.1(l are 

X = A 2 + X 

y = A r + !^ A 5 + y3A 3 + yiA (B-l) 
and the calculation of flat times l|1.14|l gives 

£ I = — Yi Xn — X™ — Y-\ Xn 

4 64 

is = §x 3 - y 3 x + ^Yi 

* - V+ 2 V ( B - 2 ) 

t 7 = Q 
2 
9 

Again we see, that (|B.2|) can be easily solved w.r.t. Yj, but the only coefficient Xo now satisfies 

^=- 3 i^-ft 5 X" - 3 -t 3 X (B.3) 

where we put t§ = for the coefficient at the "boundary" operator |55| . 
The one-point functions l|1.14(l are given for the (2, 7) model by 

|| = ™*o 8 - J^X 2 + iy 3 A 3 = -^X 5 - ~X 3 t 5 - ^X 2 i 3 

— = --yxX 3 - —X* + -Y 3 X* = -—XI - —Xth ~ ^t 3 (B.4) 

df 5 - ~512 ~ 6A YlX ° + 128 3 * ~ ~256 ~ 64 ^ ~ 128 °* 3 
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In the r.h.s.'s of i|B.4|) we already substituted the expressions for Yj in terms of times (|B.2|) . and the rest is to 
solve (|B.3|) by expanding in £3 and £5 and substitute result into (|B.4|1 . 

The scaling anzatz (|3.5f) . 13.2211 now reads 



t 2 ' ^3/2 I 

, 5 V / 

9J 7 _ ,5/2 f (l) _ X _ l/2 f (11) 

at, - h dt\ - 2 ~ * 5 ' 



(B.5) 



where we have introduced shorten notation for the derivatives over the first argument of f(ti,t2), and string 
equation (|B.3|1 turns into 

ti + ^u 4 + yu 2 + 3t 2 u = (B.6) 

for u = f( u ). Solution of JB~6|) (and then of (|B~4|) ) gives for the coefficients of f (in "broken phase" with 
nonvanishing fn), the following numbers 



= _ 7T7:V— 42, fn2 = -, fi22 = 777 V— 42, f 2 22 = — — 
90 5 35 35 

7 , 14 , 1 ^— , 6 , 9 



(B.7) 



fllll — ~ -1 „ nn V-42, fui2 — 7777' fll22 — — V— 42, fi222 — ^777' f2222 — —737V— 42 

lo20 225 50 25 175 

343 r-Jo C 196 r 49 ^fo r 28 

Tl 1 1 1 1 — V — 4z. Tl 1 1 1 9 — , tl 1 1 99 — V — 4z, Tl 1 999 — , ... 

81000 ' 3375' 2700 ' 125' 

and one can easily compute the corresponding invariant ratios, for example 

flUlfgj _ 5 f2222fll , fllllfll _ _ f2222f22 _ 1 (Ti Q\ 

f2 _ ^' f2 — _ ' f2 — _ ' f2 — ~~ o> ■■' V-°-°J 
1 112 ° 1 122 'ill '222 ° 

Ising model. Now let us turn to the "two-matrix model" Ising point (p,q) — (3,4). Here the calculation of 
residues (|1.14fl (upon vanishing t% = and £4 = 0) gives rise to 

Y 2 = \X X + 1*5 

Yi = ^X (B.9) 
2„ 2 10 



while X and Xi satisfy 



Y Q = -Xt + —X^ 
9 1 9 



*^ 2 3 ^ 2 

£l = Xr, ~\- X-t ~\- —t^X-\ 

3 27 1 9 1 
2 5 

*2 — —7:^0^1 — 2*5-^0 



(B.10) 



Upon solving the second equation for Xq , this system turns into an algebraic Kazakov-Boulatov equation |56) 
for a single function X\ 

t 1= + L X \ + -t 5 X? (B.ll) 

(2Xi + 5t 5 ) 2 27 1 9 ° 1 V ' 

which contains all information about singularities of spherical partition function for arbitrary values of magnetic 
field ti and fermion mass t§ [5*7|. 
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The one-point functions l|1.14fl are here given by 



dti 27 1 81 



. — f^X^ — —X-\Xn — —t^X^ 

, f 1 1 oi ° 1 n 1 o°U 
— — = Xf Xn ~\- t^X? 1 Xn — Xn 

dt 2 27 1 27 1 27 



(B.12) 



dU 243 1 2187 



X^\ -\- t^X^ Xt) — t^X^ — -^"n 

1 81 1 729 1 27 



Again, using 1(1.15(1 . or differentiating the first equation of <|B.12|) and using l|B.10|) one gets 

d 2 T Xi d 2 T 2X Q 



dt x dt 2 



(B.13) 



The unitary scaling anzatz 1(3.4(1 suggests that the rescaled function / depends only upon dimensionless quantities 
t ct = t 2 /t\ /6 and r e = t 5 /t{ /3 : 

T = tl /3 f( Tai T e ) (B.14) 
where two arguments of / naively correspond to the " Liouville-dressed" Ising operators a and e respectively. 



Note also, that for t 2 = the second equation of IB. 10(1 has the only reasonable solution Xo = 0, while the 
first one turns into 

t^^Xf+ 5 -t 5 X! (B.15) 

which coincides with the perturbation of the Yang-Lee (2, 5) model by quadratic term (see ((3.16(1 : note also that 
one has to identify -y* from ((B.15() with ^ from ((3.16|) . cf. with the formulas l|3.18|l and (|B.13f) \ For t 2 = 
the function l(B.14(l becomes a function of only the second argument 



T = t[ /3 f 



t 



t5_ 

1/3 



(B.16) 



Using (|B.15|I one gets for ((B.16(l 



fn = — 2 1/3 f - -— f - — 2 2/3 f - _1^2 1/3 f 



3125 
20736" 



2 2/3 



(B.17) 



while the coefficient / eeee remains undetermined from (|B.15|) . To find it one can use the first line of l|B.12|) (at 
X = 0), which gives 

/«« = 6 ~§ (B.18) 
96 

However, under reparameterization in the space of couplings (see e.g. 58 ) 

5 

Xx 



X X - -h 
125 



(B.19) 



ti 



-r' 



216 



the reduced string equation (|B.15(1 acquires the form of (analytically continued) string equation (|3.16() for the 
Yang-Lee model. Therefore, one can use further the scaling anzatz (13.19(1 with £3 = tj h of the (2, 5) critical 
point substituted by the square of the t 5 = tg Mng of the (reduced) Ising model, so that the expansion ((3.20|) 
would give the (e 2 ") correlators of the gravitationally dressed (3,4) Ising model. 
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C Residues and equations for the logarithmic potential 

Let us present here some useful explicit formulas for the two-matrix model with logarithmic potential (|3.34|1 , in 
addition to sect. 13.41 

The generating differential zdz l|3.37[) possesses four nontrivial residues at the points 

u v 

fc = *L=o = «- 7 = -- 

r us 

92 = Z \w=Vs= V + - + — ~2 = a (C.l) 

°°- = *L=oo 

°°+ = z L=s 

(one finds for the parameter (13.40ft that A = q 2 — q\), which are 



1 _ . 1 2 ru 



res„, zdz = reSoo zdz = r 



a 2 



2m' 2ni s 

11 _ u 2 ru 

— ; TeS„ 2 ZaZ — — ; TCS^, ZUZ — — - — 

2m 2m (1 — s*y s z 



(C.2) 



Taking derivative of the generating differential l|3.37|l w.r.t. the period S l|3.45|l at fixed to (|3.42[1 and introduced 
by (|3.43|l extra parameter 8, one gets the canonical holomorphic differential on l|3.36[) . 



provided by normalization condition 



d _ , dh dz 

(kj= ds zdz= -ds7m (c - 3) 

dh f dz 



Similarly, 



dh f dz [ I 1 P 



az 



dt j A VR = f A dz {w^)7n + 7n) (c - 5) 

and this completes the set of conditions, fixing all parameters in the logarithmic model of sect. 13.41 

In addition to the main cubic equation 13.49J1 . one can write a similar equation for the square of conformal 
radius (|3~1S|) 

to f X , t / A 2 



r 2 = - 



4 \t /A 2 X 

which satisfies 



(C.6) 



3 1 /A 2 7t \ 2 t /A 2 3t Uv t 2 \ t 2 
P -2\- + -- V ) p + j{- + - + 2A^-4^- 1 ') P+ 8^ it0 -'' ) = ° 



(C.7) 

The critical points of this cubic curve are at to = 0, p = 0; to = \- (1 ± 3^/1 — g^J- In the same way it is 
easy to show that the squared length of the cut 

Aru - /^+-M 2 _JL + + x-JLL (c.8) 



A 2 2A A 4 V 2 A 2 A 2 / 2A 4 X 
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also satisfies the cubic equation 

_ 3 ({t -2v) 2 M 1\„2 , *o ( l , 7i o 4 



A 4 A 2 4/ A 2 \2 2A 2 A 

t 2 (I 3 t 15t§ 4^ 2 



A 4 V 4 2 A 2 4A 4 A 4 ' A 



^(to-f)(*o-2i/))S- 
--(3to-2^) + -^(to-2^) 3 ) =0 



(C.9) 
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